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SUMMARY  PAGE 


THE  PROBLEM 

Almost  all  data  acquired  via  transducers  during  impact  experiments  must 
be  conditioned  via  some  sort  of  filter  in  order  to  optimize  signal-to-noise 
ratio  and  avoid  aliasing  of  noise  during  digitization  procedures.  To  this 
end,  rapid  design  techniques  were  developed  which  meet  nearly  all  filtering 
requirements  of  the  experimental  impact  programs  of  the  Naval  Biodynamics 
Laboratory. 

FINDINGS 

This  study  documents  the  evolution  of  a  variety  of  digital  and  analog 
filters  of  all  types  developed  over  several  years  and  presents  them  in  a 
"design  manual  format."  It  has  been  found  that  included  designs  are  adequate 
to  meet  almost  all  filtering  needs  for  data  acquired  during  biodynamic  impact 
experiments . 

RECOMMENDATIONS 

The  filter  designs  included  herein  are  all  laboratory  proven  and  can 
be  recommended  to  all  users  with  similar  needs. 
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SOME  USEFUL  FILTER  FORMS 


INTRODUCTION 

Filtering  is  an  integral  part  of  almost  all  data  acquisition  and/or  signal  conditioning 
functions.  The  most  commonly  used  filters  can  be  categorized  into  one  of  the. following 
general  types: 

1 .  Low  Pass  -  Removes  frequencies  above  some  desired  cutoff  point. 

2.  High  Pass  -  Remo\es  frequencies  below  a  certain  desired  point;  a  common  special 
case  of  this  type  is  the  "washout"  filter  which  usually  has  a  very  low  cutoff  frequency, 
and  is  primarily  used  to  remove  d.c.  bias  from  data. 

3.  Bandpass  -  A  combination  of  the  above  two  types,  it  bounds  the  signal  at  both 
ends  of  the  spectrum. 

4.  Notch  -  A  very  narrow  band  filter,  generally  designed  to  remove  a  single  unwanted 
frequency  from  the  signal. 

5.  Inverse  filter  -  Occasionally  data  is  acquired  that  was  "overfiltered"  for  some  reason, 
either  intentionally  or  due  to  bandwidth  limitation  of  the  particular  sensor  used.  Provided  the 
filtering  was  not  too  severe  and  the  filtering  transfer  function  is  known,  the  data  can  be  re¬ 
covered  by  means  of  an  inverse  filter;  such  a  device  cancels  the  effect  of  the  original  filter 
and  moves  the  cutoff  up  in  the  frequency  scale. 

There  are  vast  numbers  of  implementations  and  filtering  philosophies  available  to  the 
designer;  one  can  easily  get  bogged  down  searching  the  voluminous  literature  for  the  "opti¬ 
mum"  filter  for  a  process  and  lose  sight  of  the  fact  that  the  designer's  job  is  to  produce  an 
adequate  implementation  in  a  timely  fashion.  The  purpose  of  this  report  is  to  document  and 
provide  simple  design  information  on  a  set  of  filters  which  have  been  found  to  satisfy  most 
requirements.  Both  analog  and  digital  filters  are  discussed,  and  all  are  derived  from  the 
same  basic  second  order  transfer  function.  This  was  selected  as  the  root  function,  due  to  the 
fact  that  generally  if  a  signal  needs  filtering  at  all,  it  probably  needs  at  least  a  second  order 
(40  db/decade  slope)  filter.  Higher  (even)  order  implementations  are  possible  by  simply  cas¬ 
cading  filters,  or  in  the  digital  case,  recirculating  the  data  through  the  same  algorithm. 

ANALOG  FILTERS 

All  analog  filters  considered  herein  are  active  filters  utilizing  inverting  operational 
amplifiers.  The  latter  are  justifiable  in  that  current  technology  produces  high  quality  de¬ 
vices  at  low  cost  and  high  density  of  as  many  as  four  amplifiers  per  14-pin  DIP  package. 

The  advantages  of  using  this  active  technique  are: 

1 .  No  insertion  loss 

2.  Realization  of  complex  -  pole  2nd  order  filters  without  inductive  components 


L 


3.  Design  flexibility 

4.  Inherent  stability  of  negative  feedback  circuitry 

The  design  of  such  filters  make  use  of  the  fact  that  the  transfer  function  of  an  inverting 
operational  amplifier  is  given  simply  by  the  ratio  of  the  feedback  to  the  input  impedances, 
and  that  the  summing  junction  current  is  neglihible,  i.e.  this  junction  can  be  considered  a 
virtual  ground.  Most  analog  filters  considerea  herein  are  developed  from  a  very  useful 
configuration,  the  five  impedance  network  shown  in  figure  1;  if  one  uses  Kirchhoff's 
current  law  and  writes  nodal  equations  at  points  i  and  2  respectively,  one  obtains: 


Solving  the  equation  (2)  for  and  substituting  into  the  (1)  yields,  after  rearranging: 


This  equation  is  valid  for  any  realizable  impedance  Z,  and  is  the  basis  for  most  filters  to 
be  developed  in  subsequent  sections.  Many  other  mechanizations  are  possible,  but  the  ones 
derived  herein  have  been  found  very  useful  ar.d  simple  to  implement. 

LOW  PASS  FILTERS 

The  five  impedances  of  figure  I  can  be  defined  as  follows: 
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This  yields  the  circuit  of  figured  insertion  into  equation  (3)  and  simplification  yields: 
(Note:  Throughout  this  discussion  the  arguments  are  omitted  from  the  equations  to 
simplify  them.) 


e 

o 

e. 

i 


1  +  3  RC2 


-1 

- 5 - j 

S  +  IT  C}C2  Sz 


(4) 


where  S  is  the  Laplace  operator.  This  is  recognizable  as  the  standard  minimum  phase 
2nd  order  transfer  function  of  the  form: 


e 

o 


e. 


i 


(5) 


whereut  =  natural  (break)  frequency  in  radians/sec 
n 

^  =  damping  ratio  and  is  a  measure  of  maximum  overshoot  atto  .  By  matching 

terms  in  these  equations,  the  following  is  obtained: 


Low  pass  filters  are  usually  specified  by  the  break  frequency  and  attenuation  (in  decibels) 
at  that  frequency.  The  latter  is  related  to  the  damping  ratio  by  the  function 

An.  (db)  =  20log  ~  (8) 

obtainable  by  computing  the  magnitude  of  the  normalized  form  of  equation  (5)  atto  =  UJ>n. 

Another  parameter  of  interest  is  the  phase  shift  within  the  passband  of  the  filtered  signal. 
For  the  transfer  function  (5),  the  phase  shift  (<p)  is  given  by: 


<{>  (radians)  =  tan 


-1 


2<WUJn 

u>  z-iuz 


(9) 
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wherecuis  the  driving  frequency.  If  we  redefine  this  as  Ku^,  where  K  is  any  positive 
number,  equation  (9)  becomes: 


tan 


2  K< 


1  -  K 


1 


(10) 


For  typical  filter  damping  ratios  =  .6  to  .8)  and  the  fractions  of  the  passband 
below  approximately  20  -  30%  (Ks  .2  or  .3)  equation  10  can  be  approximated  by: 

l-K2 

since  for  small  angle$tan(p=({).  Expression  (11)  is  relevant  because  it  implies  that  for  a 
significant  portion  of  the  passband,  the  phase  shift  varies  linearly  with  frequency,  thus 
the  phase  distortion  introduced  in  this  region  consists  merely  of  a  pure  time  delay.  The 
delay  is  defined  by: 


Tj  (seconds)  =  <P  (12) 

CJU 

introducing  equation  (11)  for<P  and  again  replacing  ujwith  Ka»  ,  equation  (12)  becomes: 


2 

It  is  obvious  that  as  stated,  for  small  K,  T  ,  is  constant,  since  l-K  **  1  and: 

a  n 


Td  *  JL _  04) 

UJ 

n 

Equations  (6),  (7),  and  (8)  can  now  be  used  to  develop  an  extremely  simple  "pencil  and 
paper"  design  procedure.  Rewriting  (6)  and  (7)  in  terms  of  the  circuit  components  yields: 
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Generally,  input  design  parameters  are  break  frequency  and  attenuation  at  that 
frequency;  given  the  latter,  ^  is  computed  from  equation  (8)  .  Most  commonly  used  values 
are: 


Att.  at  break  ^ 

-  1  db  .56 

-  3  db  .707 

-  6  db  1 


In  summary,  a  typical  stepwise  design  would  be: 

Given  u)  (in  rad/sec) 
n 

a.  Select  ^ 

b.  Select  C 2 

c.  Compute  R  using  equation  (15) 

d.  Compute  C^  using  equation  (16) 

If  the  cutoff  frequency  is  specified  in  hertz  (f  )  rather  than  radians/sec,  one  replaces 
U)  with  2n  f  in  computing  R.  If  phase  shift  and/or  time  delay  are  of  interest,  they  can 
also  be  computed  using  the  given  equations.  It  should  be  noted  that  the  damping  ratio  is 
a  function  of  Cj  and  C2  only,  thus  the  cutoff  frequency  can  be  changed  via  R  without 
altering  other  filter  parameters. 

HIGH  PASS  FILTER 

The  five  impedances  of  figure  1  can  be  defined  as  follows: 


VRi 


Z  .  =00  (open  circuit) 

r» 


Z  = 


5  R2  C2  S  +  1 


(parallel  R  C  network) 


This  yields  the  configuration  of  figure  3,  and  insertion  into  equation  (3)  results  in: 


e. 

1  n 


2  2 

_ -  Rt  C1  5 _ 

2  R1  R2  cj C2  S2  +  (2R1  +  R2C2)  S  +  l 


07) 
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which  is  a  factorable  function  resulting  in  a  pair  of  real  poles: 


o  = 


2  2 

-  Rj  R2 


(1  +2  R]C1  S)  (I  +  R2C2  S) 


08) 


The  numerator  term  represents  pure  double  differentiation,  thus  attenuation  at  zero  fre¬ 
quency  is  infinite.  The  gain  increases  at  a  slope  of  +  40  db/decade  until  the  two  poles 
take  effect.  These  normally  would  be  equal  in  order  to  compensate  for  the  numerator 
dynamics;  the  first  constraint  becoming: 


2  R  C  =  R„ 

1  1  i 


C2  - 


1 


09) 


where  U^,  is  again  the  break  frequency  in  rad/sec.  An  additional  constraint  results  from 
the  determination  of  the  gain  within  the  circuit  passband.  This  is  normally  determined  by 
taking  the  limit  of  equation  (18)  as  S  approaches  infinity,  but  in  this  case  this  results  in 
an  indeterminate  form.  However,  applying  [.'Hospital's  rule  twice  gives  the  desired  para¬ 
meter; 


Passband  Gain  =  G 

P 


(20) 


combining  (19)  and  (20)  yields  another  equation  for  G  in  terms  of  the  resistors,  and  both 
of  these  expressions  must  hold:  ^ 

^  R2  (21) 

V  -np- 

A  design  procedure  can  now  be  summarized  as  before: 

Given  U)  (in  rad/sec) 
n 

a.  select  desired  G 

P 

b.  select  one  component  (usually  or  C^) 

c.  compute  other  C  using  equation  (20) 

d.  compute  R^  and  R^  from  equation  (19) 

e.  equation  (21)  can  be  used  to  verify  computation 
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This  circuit  consists  of  two  real  poles,  hence  the  attenuation  at  the  break  will  be  -6  db 
from  the  passband  gain.  If  lesser  values  are  required,  the  following  procedure  can  be 
used:  In  computing  circuit  components  R^  and  R^,  replace  ujn  byuT  the  latter  given  by 


Desired  attenuation 

-  1  db  1/3UJ 

n 

-2  db  1/2 cu 

n 

-3  db  2/3 Ut 

n 


Tbe  denominator  of  equation  (17)  is  identical  to  a  low  pass  filter,  hence  its  phase 
lag  is  obtainable  from  equation  (9)  or  (10).  To  this  one  must  add  the  numerator  contri¬ 
bution  which  is  a  constant  180  of  lead,  since  the  numerator  is  a  pure  double  differentia¬ 
tor.  Thus  the  total  phase  shift  for  the  high  pass  filter  is 

<^p  =  180°-0>Lp  (22) 

where  is  given  by  (9)  and  (10)  as  stated,  with  ^  =  1  (real  poles)  and  uu^  given  by 
equation  (19). 

NOTCH  FILTER 

A  useful  variant  of  the  five  impedance  network  is  obtained  by  defining  the  components 
as  follows: 


Z  = _ - _  •  Z  =  _ _ _ 

1  C1  S  '  L2  C2  S 

Z3  =  Z4  =  R 


Z5  C.  S 

o 

this  yields  the  network  of  figure  4;  manipulation^ similar  to  those  previously  made  yield 
the  transfer  function: 


R  C1  S 

1  +  2RC3S+R2C3  (C]  +  C2)  S2 


(23) 


this,  as  it  stands  is  of  limited  usefulness;  by  adding  another  operational  amplifier  to  the 
circuit  and  summing  the  input  to  e2,  one  develops  the  expression: 
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(24) 


which  in  terms  of  equation  (23)  becomes 

e  1+(2RC-RC.)S+R2  C  (C  +  CJ  S2 

o  _  0  1  o  I  l 

ein  1  +  2  R  C3  S  +  R2  C3  (C  +  C2)  S2 
This  equation  which  yields  the  circuit  of  figure  5  is  obviously  of  the  form 


/  S  +  S  /  2 

/  a) 

n  /  n 


i  ♦  ZW.. 5 


where  ^  ^  as  required  for  a  notch  filter.  One  might  as  well  specify^  to  be  zero, 

which  would  yield  infinite  attenuation  at(ju  .  This  is  not  realizable  in  practice,  due  to 
component  imperfection,  but  attenuation  of*1-  40  db  is  readily  achievable  with  this  specifi¬ 
cation.  This  value  of  ^  yields  the  requirement 

2  R  C3  -  R  C]  =0  (27) 

Matching  the  coefficients  of  equations  (25)  and  (26),  one  also  obtains 


W9 


=  2<j 

"  "  <29> 

The  above  three  equations  can  be  rearranged  to  yield  the  required  component  specifications 

C3  '  C1  *  2  ) 


Wn  (;  d  (C1  +  C2} 


8 


1 


A  positive  value  for  C2  requires  that^  <  .707,  an  acceptable  restriction  for  a 
notch  filter;  we  now  have  sufficient  data  to  establish  a  design  procedure: 

Given  u>n  (rad/sec)=  notch  frequency 

a.  Select  ^  (usually  (approx.)  0.6) 


b.  Select  one  component,  say  C 


c.  Remaining  components  are  computed  using  equation  (30) 
Equations  (28)  and  (29)  can  also  be  combined  to  yield  an  expression  for(^ 


^  =  C. 


V  S  (C', + 


(31) 


which  shows  that  the  damping  ratio  is  independent  of  R;  thus  R  alone  need  be  varied  to  change 
the  notch  frequency  of  the  circuit.  If  one  modifies  the  circuit  of  figure  5  to  that  of  figure  6, 
the  notch  frequency  can  be  tuned  over  a  given  range  via  the  ganged  potentiometer,  without 
affecting  any  other  filter  characteristic. 

The  denominator  phase  shift  can  again  be  obtained  from  equation  (9).  The  numerator 
contribution  consists  of  a  sudden  change  from  0°  to  +  180°  at  the  notch  frequency.  While 
this  is  not  self-evident  from  equation  (9)  for  ^  =0,  one  can  make  ^  arbitrarily  small  but 

finite,  and  see  that  the  above  holds. 

BANDPASS  FILTERS 

No  special  circuitry  is  required  to  mechanize  this  configuration.  Once  the  low  and  high 
cutoff  frequencies  are  defined,  such  a  filter  can  be  realized  by  cascading  a  high  and  a  low  pass 
filter. 

INVERSE  FILTER 

Inverse  filtering  or  "cancellation  compensation"  is  the  least  common  "on-line"  analog 
technique  since  its  function  is  to  compensate  for  an  undesirable  characteristic  of  some  signal 
path  and  presumably  a  well  designed  system  would  not  have  such  problems.  However,  there 
are  occasional  uses  for  such  a  device,  for  example  if  a  ser.sor  of  insufficient  bandwidth  is 
used  due  to  its  being  the  only  one  available.  The  required  transfer  function  is  similar  in 
form  to  that  of  a  notch  filter  except  that  unrestricted  freedom  to  choose  both  numerator  and 
denominator  parameters  are  needed.  Such  a  network  cannot  be  mechanized  via  the  five 
impedance  network.  One  convenient  configuration  can  be  obtained  via  the  "bridged-T" 
network  shown  in  figure  7.  The  transfer  impedance  can  be  readily  obtained  by  short-cir¬ 
cuiting  the  outpub 
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Z  total  = 


Z3  (Z4  +  Z5} 
Z3+Z4  +  Z5 


Z 


1 


+ 


+ 


Z3  <Z4  +  V 
Z3  +  Z4  +  Z5 


(32) 


Z  total  = 


Z1  Z2  Z3  *  Z1  Z2  Z4  +  Z1  Z2  Z5  +  Z1  Z3  Z4  +  Z|  Z3  Z5 _ 

Z1  Z3  +  Z1  Z4  +  Z1  Z5  +  Z2  Z3  +  Z2  Z4  +  Z2  Z5  +  Z3  Z4  +  Z3  Z5 


(33) 


If  one  defines 


z  =  Z  =  — - _ 

2  j  C,  S 


Z4  C„  S 


VR1 

The  network  of  figure  8  is  obtained  and  the  transfer  impedance  becomes,  after  some  tedious 
but  straightforward  manipulations: 


Z  total  =  R„ 


1  +  A  S 


2  2 

B  S  +CS+1 


(34) 


where  A  = 


B  = 


C  = 


2  R1  C]  C2 

2  Cj  +  C2 

R,  R2  C,2C2 
2C,  +C2 

C1  (2Rt  C2  +R2  Cp 


y 


(35) 
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Since  the  overall  transfer  function  of  an  operational  amplifier  is  given  by  the  ratio  of  the 
feedback  to  the  input  impedances,  the  insertion  of  identical  networks  for  the  two  impedances, 
as  shown  in  figure  9,  results  in  the  form: 


(1  +  a's)  (BS2  +DS  +1)  (36) 

(1  +AS)(b'  S2  +d‘  S  +  !) 


it  is  evident  that  if  A  is  constrained  to  be  equal  to  A  ,  a  transfer  function  of  the  desired  form 
obtains. 

In  inverse  filtering  the  numerator  coefficients  B  and  C  are  forced  to  match  those  of 
the  undesired  filter  function,  whereas  B  and  D  are  made  to  be  a  multiple  of  B  and  D,  hence 
extending  the  overall  system  bandwidth.  Thus  one  can  express  (36)  in  the  form 


e 

o 

e. 
i  n 


iL  +  2-£l  +  1 

<  -n2 


S2  + 2  +  1 

K  ujn  Ku^ 


(37) 


where  K  is  a  positive  number  relating  the  new  break  frequency  to  the  one  being  cancelled. 

Since  the  damping  ratios  may  or  may  not  be  the  same,  they  are  maintained  as  distinct  quantities. 
Invocation  of  equations  (35),  (36),  and  (37)  is  sufficient  to  define  the  circuit  components  in  terms 
of  the  filter  parameters,  but  these  equations  are  unwieldy.  Great  simplification  is  possible  by 
introducing  a  couple  of  constraints 


A'  =A  =  D/10 


R 


2 


I 


=  R 


2 


The  first  of  these  moves  the  first  order  terms  to  be  cancelled  to  a  high  frequency  region,  where 
possible  mismatches  would  be  inconsequential;  it  does  limit  K,  but  values  up  to  6  are  allowed, 
thus  covering  the  range  of  practical  interest.  The  second  constraint  merely  assures  unity  d.c. 
gain. 


Implementation  of  the  above,  coupled  with  insertion  of  the  basic  parameters,  yield  the 
component  specifications 


n 


R 


1 


C 


1 


I 


10 _ 

1  n 


C 


2 


R 

2 


( _ s - 

^Ku^OO^-K^) 


R  is  unrestricted  and  user  selectable.  Greater  computational  simplicity  results  from  the 
inverse Torms  for  Rj  and  R^ 


There  are  certain  conditions  implied  by  the  R  .  and  C 2  specifications  since  these  components 
must  be  positive.  The  most  restrictive  is  that  For  C  2  where  it  is  required  that 

I0C2>  KC1 

which  at  worst  (for  K=6),  requires  that  ^  ,  a  condition  that  is  easily  met  in  almost  any 

practical  situation.  A  design  procedure  for  an  inverse  filter  is  now  possible 

Given 

UUn  (rad/sec),  K,  (K  4  6)  ^,^2 

I 

a.  Select  R2  (&2  =  ^ 

b.  Compute  remaining  components  using  equation(s)  (38). 
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It  should  be  noted  that  for  low  frequencies  (u>  <C  10),  either  1^  or  C ^  must  be  quite 
large.  This  is  unavoidable  in  this  network  configuration,  thus  if  a  very  low  frequency  de¬ 
sign  is  required,  alternate  forms  (ref.  1)  may  be  more  appropriate  if  these  large  component 
values  are  unacceptable. 

The  network  phase  shift  is  given  by 
♦total  =  V^2 

whereqv  and  <|>«  are  obtained  by  applying  equation  (9)  to  the  numerator  (<10  and  denomin¬ 
ator  (<t>2)  of  equation  (37). 

DIGITAL  FILTERS 


Frequently  data  that  has  been  acquired  and  transformed  into  digital  format  is  found 
to  require  additional  filtering  of  any  of  the  types  discussed  in  the  Introduction.  In  the 
analog  domain,  each  filter  type  was  developed  using  special  dedicated  circuit  in  order  to 
obtain  an  economical,  yet  simple  practical  design.  Since  a  digital  filter  is  generally  a 
software  implementation  of  an  algorithm,  all  five  filter  types  can  be  realized  from  the 
same  basic  mechanization  of  the  most  general  form  of  the  second  order  equation,  depicted 
by 


S  +  1 


S  +  1 


(41) 


Equation  (41)  can  be  modeled  via  a  "simulation"  block  diagram  as  if  it  were  to  be 
programmed  into  an  analog  computer.  More  than  one  configuration  is  possible,  but  herein 
the  "M"  method  (ref.  2)  will  be  used.  This  technique  results  in  an  analog  computer  block 
diagram  as  shown  in  figure  10;  in  this  most  general  representation,  an  inverse  filter  is  des¬ 
cribed,  but  a  notch  is  readily  achieved  by  simply  inserting  the  appropriate  coefficient  values 
(i.e.^  j  =  0,  lUj  =  u>2  ).  For  a  low  pass  filter  the  numerator  of  equation  (41)  needs  to  be 
unity,  which  is  obtainable  if  we  let CU^  approach  infinity.  If  this  is  done,  it  is  evident  from 
figure  10  that  EQ  =  Xj ;  thus  a  low  pass  filter  is  also  obtained  from  the  mechanization  of 
figure  10,  by  using  Xj  as  the  output. 

While  less  obvious,  the  high  pass  configuration  is  also  obtainable  from  this  diagram. 
Such  a  filter,  with  unity  gain  in  the  passband  is  described  by 


(42) 


- 


From  figure  10  one  oan  work  backwards  from  Xj  to  the  point  labeled  X^ ,  and  observe  that 

.2 


X3  = 


1 


X1  " 


cu. 


‘l 


(43) 


Ui2  w2 

But  X.  has  already  been  shown  to  be  the  low  pass  output,  thus  equation  (43)  is  equivalent 
to  (42)  and  one  obtains  the  high  pass  configuration  by  defining  X^  as  the  output. 

Algorithm  Development.  The  derivation  of  a  digital  filter  can  be  obtained  via  the  state 
space  equations  that  describes  the  system  shown  in  figure  9.  The  state  equations  are  discret¬ 
ized  and  then  solved,  and  the  solution  iterated  over  time.  This  technique  was  originally 
described  in  reference  3,  and  is  briefly  repeated  herein  for  completeness.  From  the  figure, 
the  vector-matrix  filter  equations  are  readily  obtained 


X, 


L  2  J 


ix>2 ' 


1 


2<2  “Y 


M 

'  0 

1 

X 

ro 

I _ 

+ 

2 

0)2 

—  m 

in 


E  ^  is  obtainable  as  a  linear  combination  of  X,  and  X,. 
out  1  2 


(44) 


ai. 


"out 


Ui 


E  (1  -u>22)  x 

in  + - j -  I  + 


W 


1 


-2C2U,2)  X 


U>, 


2  (45) 


The  above  vector-matrix  differential  equation  is  of  the  general  form 
X  =  Ax  +  Bu 
whose  solution  is  given  by: 


(46) 


X(t)  =  eAt  X(o)  +  e  _ArBu  (r  )  dr  (47) 

Since  the  data  of  interest  is  discrete  (due  to  the  digitization  process)  the  solution  for  the 
discrete  form  of  subject  set  of  differential  equation  is  needed.  A  recursive  relation  ex¬ 
ists  that  allows  the  computation  of  the  system  state  at  time  (K  +  1)T  based  on  knowledge 
of  the  state  at  time  KT;  K  is  a  running  index  and  T  is  the  sampling  interval.  The  recursive 
relation  is  given  by  a  formula.  (Reference  4) 

X(K  +  1)  T  =  G(T)  X  (KT)  +  H(T)  u  (KT)  (48) 
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where  G(T)  =  eA*;  H(T)  =  f  ^  e^T  dT  B  (49) 

«/  o 

The  solution  of  the  subject  difference  equation  reduces  essentially  to  computing 
e  and  its  integral  over  T,  a  tedious,  but  straightforward  task  which  is  rendered  simpler 
by  the  following  substitutions: 


(;  cu  =  a.  cu=cu  Ji  - 

?  •  2  s  i 


*2  2  '  2 
The  resulting  solution  then  becomes 


(50) 


X,  (K+l) 

= 

A  A 

*11  *12 

’x  1  (K) 

+ 

CU  ^  F 

2  hl2 

x2  (K+l) 

i 

<N 

< 

CN 

x2  (K) 

i - 

J 

n: 

*n 

N> 

ro 

1 - 

E.  (K)  (51) 

in 


where 


=  e  aT(cosu>T  +  sin  CUT) 


(52) 


A12=  iS  (e  aTs5n«»T) 


(53) 


A21 


(°2^U,2L[e-°Ti,nU)T) 


A 22  ~  ®  OT|cos  CUT 


jcos  CUT  -  aXu^in  CUt| 


(54) 

(55) 


-aT 


F«2  2  2 

u  cu(cr  +  cu  ) 


(  asin  CUT  +  CUcos  CUT)  +  — s - 9 

,  fll* 


a*  +  cu 


(56) 


F22  A12 


From  the  original  equation,  (K)  is  inferred  to  be 


E  (K)  = 

out 


CU' 


cu 


(-=r) 


(57) 


T-Ein  (K)+  ,-TT  X!(K)  + 


(g^cu,  -2^2cu2) 


CUi 


x2(K) 


(58) 
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Therefore,  the  digital  algorithm  consists  of  mechanizing  the  three  equations 

for  X  (K  +  1),  X-  (K  +  1)  and  E  (K).  Required  inputs  are  the  sequence  to  be 

filtered  E.^  (K),  and  the  initial  s?ates  (0)  and  X^  (0)  if  known.  If  the  input 

sequence  Vo  be  filtered  is  non-zero  at  the  start  of  the  process,  it  is  usually  sufficient 

to  set  X  (0)  =  E.  (0). 

I  in 

The  substitutions  (50)  were  made  for  convenience  in  computing  certain  inverse 
Laplace  transforms.  Reinserting  the  original  values  into  equations  (52)  through  (57) 
simplifies  them  considerably;  additionally,  (54)  through  (57)  can  be  expressed  in  terms 
of  (52)  and  (53)  as  follows: 


An  = 

e  ^2^*2  ^  (cosJ3u)2 

T  + 

<2 

p  sin  pa>2  T) 

(59) 

A12  = 

sin  3 

2tW2T> 

(60) 

V'  -  <2* 

A21  = 

CN 

<r 

CN 

CN 

3 

1 

(61) 

li 

CN 

CN 

< 

-2^2“2  A12  +  All 

(62) 

F>2  - 

-V  <'-Au> 

UJ 

(63) 

(64) 


Thus  the  equations  to  be  programmed  and  solved  are 

X,  (K  +1)  =  (X}  (K)  -  E.n  (K)]  An  +  A]2X2  (K)  =  EJn  (K)  (65) 

X2(K  +  1)=  (E.p(K)-X1  (K))  (w22  Al2)+(Air2C2ui2  A^)  (K)  (66) 


A  typical  design  procedure  would  be: 


Given 


U)| ,  C|  (if  applicable) 

w2,  ^2 
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a.  Select  iteration  interval  T 

b.  Compute  A^  j,  (they  are  constants  for  a  given  filter) 

c.  Insert  into  iteration  algorithm 

Different  filter  types  are  obtained  as  follows: 

Low  Pass  -  Solution  of  equations  (67)  and  (68)  is  sufficient,  since  the 
filtered  output  is  (K). 

High  Pass  -  As  previously  stated,  is  given  by  Xg(K)  which  is  from  figure  10 

2C2 

X3  (K)  =  E.n  (K)  -  Xt(K)  -  X2  (K) 


(67) 


Notch  -  Given  by  equation  (60)  simplified  by  u>^  =  uu2  and  { ^  =  0 

2  <2 


Eo<K>=Ein<K>-^-*2 


(68) 


Inverse  -  Given  by  equation  (58),  unaltered. 

Bandpass  -  As  in  the  analog  case,  realized  by  circulating  data  through 
separate  high  and  low  pass  filters. 

Iteration  Interval  and  Possible  Simplifications 


The  iteration  interval  T  has  thus  for  not  been  quantified;  it  should  however 
be  sufficient  small  to  avoid  aliasing  the  data,  including  the  noise  components  to  be  filtered. 
It  is  actually  possible  to  heuristically  quantify  T  in  terms  of  the  filter  parameters  which  leads 
to  normalization  of  the  coefficient  (A.  ^ ,  A  equations  and  opens  the  door  for  some  very 
useful  simplificationsof  the  filtering  algorithm. 


The  following  conjecture  is  proposed:  If  a  signal  needs  to  be  filtered  at  all, 
it  must  have  frequency  components  significantly  higher  than  the  filter  cutoff,  and  the 
assumption  of  frequencies  10  times U>2  is  not  unreasonable.  In  addition,  the  sampling 
frequency  must  be  high  enough  to  reasonably  duplicate  this  component.  A  minimum 
specification  then  can  be  developed  where  the  sampling  frequency  should  be  at  least 
50  times  the  filter  cutoff. 


Defining; 


CO^  =  sampling  frequency  =  -y 
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one  obtains,  as  a  result  of  the  above  conjecture,  the  relationship: 


—  <  0.02 

U)  ~ 

s 

It  should  be  noted  that  Hi  =  — -  is  not  a  sinusoidal  frequency  in  the  strict  sense,  but 
it  is,  nevertheless  convenient  to  use  this  notation  for  the  sampling  rote.  Substituting 
(69)  into  (59)  and  (60)  yields: 


(70) 


A 


11 


=  e 


(71) 


(72) 


However,  the  term  0  =  ^1  will  always  be  less  than  unity,  and  equation  (70) 

is  an  upper  limit  for  the  sine  and  cosine  values  in  equation  (71)  and  (72).  For  small 
angles  sin  X  =  X  and  cos  X  =  1 ,  and  one  can  then  simplify  (71)  and  (72). 

OJ0 

(73) 


(74) 
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These  equations  are  accurate  to  better  than  0.02%  subject  to  the  constraint  of  (70) 

Even  when  the  ratio  of  is  as  large  as  0.1,  the  equations  retain  an  error  of  <  0.5 /o, 

Equation  (73)  can  be  further  examined  in  terms  of  the  parameters  ^  and  UtyU>$ 
Confutation  of  Aj  1  in  the  ranges  ***2/^5  -  *02  an(*  0*  I  ^  "  ^.8  yie*ds  values  for  Aj  j 


ranging  from  0.99  to  0.998.  Hence,  is,  for  all  practical  purposes,  unity,  resulting 
in  the  following  simplified  set  of  filter  equations,  which  in  addition  to  the  above  simpli¬ 
fication  have  had  replaced  by  a  new  variable  =  ^  ^'s  norTnal'zes 

B^2  ar>d  006  can  pre-compute  it  as  a  function  of  the  damping  and  frequency  ratios: 


Xj  (K  +  1)  =  Xj  (K)  +  X2  (K) 


(75) 


X2  (K+l)  = 


[E;n(K)-X,(K>] 


U>2B12  +0-2C2B12)X2(K) 


-r 

***2  ^2  tu 

8  =  _ _  6  S 

B12  U)  ® 
s 


(76) 

(77) 


A  non-transcendental  approximation  for  (77)  is  given  by: 


Bl2  *0  -  .1<2> 


uj2 


u> 


(78) 


The  accuracy  of  this  linear  approximation  is  biased  in  favor  of  the  smaller  values  of  ^2» 
the  less  stable  region  for  the  algorithm. 

The  areas  of  applicability  for  the  approximate  forms  will  be  discussed  in  some  detail 
in  the  experimental  results  section.  Where  useable,  these  forms  are  obviously  simpler  to 
implement  and  the  coefficients  readily  calculated.  This  is  particularly  handy  for  "on-line" 
microcomputer  applications:  These  devices  frequently  do  not  have  a  hardware  multiply  and 
divide  capability,  and  since  one  usually  has  some  selection  flexibility  on  filter  parameters, 
it  is  often  possible  to  constrain  these  to  integral  powers  of  2  (or  summations  thereof).  This 
reduces  all  computations  to  no  more  than  add  (or  subtract)  and  shift  operations,  thus  greatly 
reducing  execution  time.  The  exact  equations  and  related  coefficients  are  summarized  in 
Table  1;  Table  2  consolidates  the  various  approximation  levels  discussed. 
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EXPERIMENTAL  RESULTS 


Analog  mechanizations  for  low  pass,  high  pass,  and  notch  filters  are  illustrated 
in  figures  11  -  13.  Figure  14  shows  a  cascaded  combination  of  a  150  Hz  low  pass 
filter  and  a  400  Hz  notch. 

Both  approximate  and  exact  digital  low  pass  and  high  pass  configurations,  are 
illustrated  in  figures  15  and  16.  Figures  17  and  18  depict  a  low  pass  configuration 
with  low  damping  ratio  =  .2),  which  while  not  a  useful  filter  form  does  indicate 
the  applicability  of  the  algorithm  to  the  more  general  problem  of  transfer  function 
simulation . 

Figures  19  -  22  illustrate  the  effect  of  the  various  approximations  and  sampling 
rates  on  the  relative  quality  of  notch  achievable.  Finally,  the  applicability  of  the 
digital  implementation  to  the  general  2nd  order/2nd  order  transfer  function  (inverse 
filter)  is  illustrated  in  figures  23  and  24. 

It  is  evident  that,  with  the  exception  of  the  notch  filter,  in  almost  all  cases  it 
makes  little  difference  which  form  of  the  equations  and  coefficients  are  used;  conse¬ 
quently,  (for  Ts  .02)  one  might  as  well  use  the  simplest  implementation  possible. 

The  notch  filter  is  not  unexpectedly,  somewhat  more  sensitive  to  approximations; 
here  the  designer  must  weigh  the  filtering  requirements  versus  design  simplicity. 

One  additional  point  should  be  made  with  respect  to  the  digital  implementations: 

In  addition  to  the  effects  of  algorithm  mechanizations,  the  designer  should  be  aware 
of  the  contribution  that  the  truncated  arithmetic  of  the  computer  itself  might  make, 
a  point  not  directly  addressed  herein,  but  implicitly  evident  in  the  notch  filter  illustra¬ 
tions. 

The  theoretical  algorithm  (equation  68)  yields  an  infinite  notch  which  was  obvi¬ 
ously  not  achieved  just  as  it  was  not  in  the  analog  case.  Whereas  the  analog  degradation 
is  attributable  to  component  imperfections,  the  digital  one  is  caused  by  computational 
"imperfections."  The  illustrated  filters  were  implemented  on  a  16  bit  minicomputer 
utilizing  FORTRAN  and  single  precision  arithmetic.  The  utilization  of  longer  word 
length  computers  and/or  multiple  precision  arithmetic  will  significantly  improve  the 
quality  of  the  notch.  The  relative  ranking  of  the  various  approximate  forms  would 
likely  remain  the  same  at  this  higher  performance  level. 
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FIGURE  2  LOW-PASS  FILTER 
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FIGURE  5  NOTCH  FILTER 
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FIGURE  9  INVERSE  FILTER 


IX1(K+1)=  {Xj( K ) -  E,n(K)}au  +  A12X2(K)+  Ein(K) 

x^k+d^IejnCkj-x^k^Ia^+Iah-z^^a^Jx^ 

(au=  e"^T{cOS(fJu>2T)4^2/3  SIN(pu>2T))} 
A12=  l/^ple^^siNO^T)]} 


p 


LOW  PASS:  Xj(K) 

HIGH  PASS:  E|N(K)~ 

NOTCH:  E|n(K)-2?2/^x2(K)} 


INVERSE^ujf/w^EJK)}-.  > i 
Uij=2npj 
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T=  l/^s  W5  =  SAMPLING  FREQUENCY 


TABLE  1  DIGITAL  FILTER  EQUATIONS 


APROXIMATE  COEFFICIENTS 
u*2/ws<0.02 

Iau= 

I 

(A12=e“<;2UJ2/UJs|i/ujs| 

EQUATIONS  SAME  AS  TABLE  1 


APPROXIMATE  EQUATIONS 


II 


X1  (K+l)=  XX(  K)  +{B12/oi2({x2(K)| 
X2(K+l)={E|N(K)-X1(K)}ui2B12  +  |l-2C2B12  }x2(K) 


|  1-0*K2  |uj2/ws 


FILTER  TERMS  SAME  AS  TABLE  1 

TABLE  2 

APPROXIMATIONS  FOR  DIGITAL  EQUATIONS 
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FIGURE  13-NOTCH 


FIG.  14-COMBINATION  LOW  PASS  A  NOTCH 


FREQUENCY  (NORMALIZED) 


FREQUENCY  (NORMALIZED) 


FREQUENCY  (NORMALIZED) 


FREQUENCY  (NORMALIZED) 


1 ,  A  ^  2  Exact 

Exact  Aj2  -  Linear  Approx. 
1  A. «  =  Linear  Approx. 


FREQUENCY  (NORMALIZED) 


FREQUENCY  (NORMALIZED) 


FREQUENCY  (NORMALIZED) 


FREQUENCY  (NORMALIZED) 


Appendix 

An  Alternate  High  Pass  Form 


The  high  pass  filter  discussed  in  the  main  body  of  the  report  yields  real  poles  only. 
At  times  one  may  need  complex  poles  and  in  this  instance  a  form  very  similar  to  the  low 
pass  five  impedance  network  can  be  obtained  by  defining: 


Z2  =  R, 


Z5  =  R2 


This  yields  the  network  of  figure  A-l  and  the  resulting  transfer  function  is  given  by 


2  2 

e  R.  R9  <T  SZ 

o  12 

—  -  - * — x - 

ein  Ri  R2  C  S  +  3R1  CS  +  1 


with  parameters 


Passband  Gain  =  1 

1 


UJ 


"  R1  *2  °2 


?  =  1 


.5^7* 


The  similarity  of  the  equations  to  those  for  the  low  pass  filter  is  obvious,  with  the 
resistors  and  capacitors  interchanging  roles.  The  component  definitions  are: 


c  = 


1.5  R, 


u> 


2  25 

R|  =  — r2 
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and  a  proposed  design  procedure 
Select  l?2 

Compute 


Compute  C 

All  phase  shift  and  time  shift  equations  developed  for  the  low  pass  case  apply 
without  change  except  that  the  phase  is  positive  (leading)  and  the  time  shift  is  an  advance. 


f  R2=SELECTED 


MEGOHMS 

AND 

MICROFARADS 


I  Rj-2-25R2/^2 
|C=  C/l.5R2tuN 


Gp=l 


=  20fN 


FIGURE  A-l  ALTERNATE  HIGH  PASS  FORM 
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